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Abstract 



The neutron is electrically neutral, but its substructure consists of charged quarks so it may 
have an internal charge distribution. In fact it is known to have a negative mean square charge 
radius (MSCR), the second moment of the radial charge density. In other words the neutron has 
a positive core and negative skin. In the first Born approximation the neutron MSCR can be 
simply related to the neutron-electron scattering length b ne . In the past this important quantity 
has been extracted from the energy dependence of the total transmission cross-section of neu- 
trons on high-Z targets, a very difficult and complicated process. A few years ago S. A. Werner 
proposed a novel approach to measuring b ne from the neutron's dynamical phase shift in a per- 



fect crystal close to the Bragg condition. We are conducting an experiment based on this method 
""^i . at the NIST neutron interferometer which may lead to a five-fold improvement in precision of 



b ne and hence the neutron MSCR. 



5 ■ 1 Introduction 

^ ' 

The neutron has zero net electric charge, but it is composed of charged quarks so it may have 
an internal charge distribution. If so the neutron will interact electrostatically with other charged 
particles such as electrons. In a scattering experiment it is useful to parameterize the neutron 
charge distribution using the Sachs electric form factor Ge(Q 2 ), the spatial Fourier transform of 
the neutron charge density in the Breit (brick wall) frame, i.e. in the frame where the neutron's 
initial momentum is Pi = —Q/2 and its final momentum is pf = +Q/2 pQ. Writing the electric 
form factor as a multipole expansion in Q 2 we have: 

G E (Q 2 )=q n -\(rl)Q 2 + ... (1) 
o 

At Q 2 = the electric form factor is simply the neutron's net charge: q n = 0. In the low Q 2 limit 
Ge(Q 2 ) is dominated by the second term which contains the neutron mean square charge radius 
(MSCR) (r 2 ), the second moment of the neutron charge density: 

(r 2 n ) = J P(r)r 2 d 3 r. (2) 

From Q we see that (r 2 ) determines the slope of the form factor at Q 2 = 0: 

dQ 



«>=-6^> , . (3) 
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In the subfield of nuclear physics concerned with understanding nucleon structure, one measures 
Ge{Q 2 ) as a function of Q 2 by scattering electrons from light nuclei at various energies and compares 
it to the predictions of theoretical models. For this program (J3J) is an important constraint on the 
shape of Ge{Q 2 )- 

Theoretically there are several ways of understanding how a nonzero (r 2 ) might arise. The 
simple SU(6) quark model predicts Ge(Q 2 ) to be identically zero, and hence (r 2 ) = 02. The 
known A— n mass splitting implies a spin-dependent quark interaction that breaks SU(6) symmetry. 
Such a force would also tend to pull the up quark to the center of the neutron and push the down 
quarks to the edge, producing a nonzero Ge(Q 2 ) and specifically a negative (r 2 ) [5J- In another 
picture every neutron spends part of its time as a virtual proton-pion pair [3]. In this state the 
lighter ir~ orbits the heavier proton, making a negative (r 2 ). Recent theoretical efforts to calculate 
(r 2 ) produce a negative sign but a large uncertainty on its magnitude jEHEHZj. 

In an experiment, the total coherent scattering length for neutron-atom scattering (far from 
any resonances and neglecting small effects such as the neutron's electric polarizability) is given 
by: 

fecoh = b N + Z [1 - /(g)] b ne (4) 

where b^ is the nuclear scattering length and b ne is the neutron-electron scattering length, respon- 
sible for the electrostatic interaction between the neutron and a charged particle {e.g. electron). 
The first term in brackets corresponds to the nuclear electric charge and the second term to the 
atomic electrons. The form factor f(q) is the spatial Fourier transform of the atomic electron 
charge distribution for momentum transfer q. 

The net charge of the particle and its internal charge distribution both contribute to b ne . 
For neutron- atom scattering at low momentum transfer the second term in Q dominates. From 
scattering theory, using the 1 st Born approximation, one obtains: 

(r 2 n ) = 3a (—) b ne = (86.34 fm)-6 ne . (5) 
\m n J 

In recent experimental determinations of b ne , total transmission cross sections are measured for 
epithermal and thermal neutrons on lead and bismuth targets [HI 03 EH- These studies, and also 
the earlier work of Krohn and Ringo jllj using noble gas targets, contribute to the Particle Data 
Group's recommended value of (r 2 ) = —0.1161 ±0.0022 fm 2 [T^]. This includes a scale factor of 1.3 
to account for a minor disagreement. There is a much larger and still unexplained disagreement 
between measurements of b ne from transmission measurements and older determinations neutron 
diffraction on crystals [13]. 

2 A New Method 

Recently S. A. Werner proposed a new approach for a precision measurement of b ne using neutron 
interferometry. The basic experimental setup is illustrated in Figure ^ which shows a standard 
LLL perfect crystal silicon neutron interferometer and the usual phase flag (e.g. fused silica). The 
neutron beam is split into two paths by Bragg diffraction, and the relative phase shift caused by 
different path lengths in the phase flag produces interference fringes in the neutron counters C2 and 
C3. Detailed introductions to the theory and operation of neutron inteferometers can be found in 
[141 1151 FTH] . A separate perfect crystal silicon sample, cut from the same ingot as the interferometer, 
is placed across both downstream paths. This sample can be rotated about all three axes. It is 
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Figure 1: A sketch of the experimental setup. 



aligned precisely to the Bragg condition (parallel to the interferometer blades) and then rotated by 
a small misset angle 59 as indicated. The neutron phase as a function of 59 over a range of ±10 
arcsec from Bragg is then measured. 

Very close to Bragg the neutron obtains a dynamical phase shift that originates from the 
transition to the a and (3 solutions for the internal wave vector. According to dynamical diffraction 
theory this is given by |17j : 



where the minus (plus) sign corresponds to the a {(3) solution. Here ujj is the probability for Bragg 
scattering per unit path length, related to the crystal structure function Fh for this reflection by 
vn = -Fff A/Vceii- The parameter y is a scaled misset angle: 



D is the crystal thickness, and 9b is the Bragg angle. Because the crystal sample extends through 
both beam paths the kinematic phase cancels and the dynamical phase adds, producing a net 
relative phase: 



with the plus (minus) sign corresponding to positive (negative) misset angles. A plot of A<I> vs. 59 
is shown in Figured 

The sample crystal is suspended on a shaft connected to a 10 cm lever, which is moved by 
a self-calibrated PZT nanopositioner that steps over a range of 25 fim with 0.1 nm precision (P- 
753, Physik Instrumente, Karlsruhe, Germany JH])- The shaft rotates on a nylon flexure bearing, 
there are no sliding or rolling contacts, so the crystal rotates smoothly and the misset angle can be 
measured to 10~ 4 arcsec. 

The neutron interferometer can measure the net relative phase shift A<i> to within 1° in a 
single interferogram (in about 1000 s) The silicon lattice constant, neutron wavelength, and crystal 
thickness D are all known at the 10 -5 precision level. Therefore by fitting the measured phase shift 
as a function of misset angle to Equation |H] one can with some effort determine ujj, and hence Fh, 




(6) 




(8) 
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Figure 2: The theoretical relative phase A<I> as a function of the misset angle 
calculated from Equation |HJ 



from Bragg, 



to an absolute precision of about 3 x 10 5 . Near the Bragg condition the structure factor is related 
to the coherent scattering length (see@J) by: 

F H = ^32 (b N + Z [1 - f{H)\ b ne ) . (9) 

In our case f(H) = 0.7526 |19j is the form factor calculated for the reciprocal lattice vector 
H corresponding to the (111) reflection in silicon. In Equation E3 6tv is about 240 times larger 
than Z [1 — f(H)] b ne . By rotating the crystal off Bragg far enough that the dynamical phase is 
negligible we can measure 6 at in situ with a relative uncertainty of less than 10~ 5 , limited only 
by the uncertainty in the phase measurement, and subtract it from Equation |UJ This will allow 
a determination of b ne and hence (r^) with a net uncertainty of about 0.7%, more than a factor 
of three improvement from the present accepted value, and using a technique systematically much 
different from other precision measurements. If the experiment with silicon is successful we may 
repeat it using a perfect germanium crystal, improving the sensitivity by the increase in Z: a 
factor of 2.3. This may lead to a 0.3% measurement of b ne and (r%), more than a factor of five 
improvement. 

We note that this method is related to the experiment of Shull which first demonstrated 
Pendellosung interference fringes in neutron dynamical diffraction [201 - ^ ne structure factor Fh 
was precisely determined from the fringe separation and 6 co h was obtained for the (111) reflection in 
silicon via Equation EO The main advantages of our approach are (1) an improved overall precision 
due to the interferometric measurement of the neutron phase and (2) the ability to make a relative 
subr action of b^ from Equation ||iti situ. 
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